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1.  INTRODUCTION 


The  calculation  of  infrared  (IR)  radiance  and  transmittance  spectra  is  an  important 
ingredient  in  many  areas  of  atmospheric  science.  These  include  modeling  the  atmospheric 
energy  budget  for  global  climate  change  predictions,  analyzing  data  from  remote  sounding 
experiments,  and  understanding  molecular  excitation  and  production  processes.  Typically, 
a  large  number  of  calculations  are  required  that  span  an  extensive  spectral/spatial  domain. 
For  the  upper  atmosphere,  non-local  thermodynamic  equilibrium  (non-LTE)  effects  must 
also  be  treated. 

Available  computer  codes  for  atmospheric  IR  radiation  that  have  gained  wide  use 
include  several  developed  and  maintained  by  the  US  Air  Force.  LOWTRAN1  and 
MODTRAN2  are  rapid,  LTE  codes  for  low-  to  moderate-resolution  spectra.  FASCODE3 
is  a  high-resolution,  line-by-line  code  that  can  be  used  in  both  LTE  and  non-LTE 
applications;  however,  the  non-LTE  molecular  state  populations  must  be  externally 
generated.  Until  recently,  a  rapid  and  unified  standard  code  for  first-principles  calculations 
of  non-LTE  atmospheric  background  radiation  has  been  lacking. 

This  report  briefly  describes  technical  aspects  of  a  new  Air  Force  non-LTE  code  for 
atmospheric  IR  spectra  that  has  sufficient  accuracy,  speed,  and  resolution  (0.5  cm'1)  for  a 
wide  range  of  applications  and  provides  a  comprehensive  treatment  of  all  important 
molecular  species.  SHARC,  the  Strategic  High-altitude  Atmospheric  Radiation  Code,4"6 
calculates  emission  and  transmittance  spectra  in  the  2-40  /xm  region  for  arbitrary 
line-of-sight  (LOS)  paths  between  50  and  300  km.  It  incorporates  the  IR  bands  of  NO, 
CO2,  O3,  H2O,  OH,  CO,  and  CH4  found  in  the  quiescent  atmosphere,  including  minor 
isotopic  bands  of  C02  and  H20.  It  also  accounts  for  auroral  production  and  excitation  of 
C02,  NO,  and  NO+  caused  by  the  flux  of  energetic  solar  electrons.  Proton  fluxes  are  not 
presently  modeled  but  could  be  accommodated  in  a  similar  manner. 

To  account  for  non-LTE  effects,  molecular  vibrational  state  populations  are 
calculated  from  first  principles.  As  in  other  non-LTE  atmospheric  models  (e.g., 
Wintersteiner  et  al.,7  Lopez-Puertas  et  al.8,9)  steady-state  kinetics  are  assumed  for  the 
quiescent  atmospheric  processes  of  collisional  excitation,  de-excitation,  energy  transfer, 
radiative  decay,  illumination  by  the  sun,  earth  and  atmosphere,  and  chemical  production. 
In  SHARC,  additional  production  and  excitation  mechanisms  resulting  from  electron 
deposition  are  included  using  a  time-dependent  kinetic  model.  Rotation -translation 
equilibrium  is  assumed. 


The  calculated  vibrational  state  populations  are  input  to  a  LOS  spectral  radiance 
module,  which  uses  a  line-by-line  (LBL)  equivalent  width  formalism.  If  desired,  the 
populations  may  be  fed  to  an  external  spectral  model  such  as  FASCODE,  which  can 
provide  better  resolution  and  accuracy,  although  at  a  considerable  time  penalty.  Our  own 
comparisons  indicate  nearly  a  two-order-of-magnitude  improvement  in  speed  using  the 
equivalent  width  approach  compared  to  the  standard  LBL  grid  method,  with  typical 
radiance  errors  of  less  than  10%. 

SHARC  is  written  in  FORTRAN-77  and  has  a  modular  structure  to  facilitate 
troubleshooting,  modifications,  and  additions  to  the  kinetic  and  spectroscopic  databases  and 
molecular  species.  A  menu-driven  user  interface  is  included.  In  addition  to  its  usefulness 
for  strategic  background  calculations,  SHARC  should  be  valuable  for  atmospheric  remote 
sounding  and  energy  transport  applications.  The  code  is  available  for  use  by  the  scientific 
community,  and  may  be  obtained  from  the  Phillips  Laboratory/GPOS  (formerly  the  Air 
Force  Geophysics  Laboratory).10 
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2.  CODE  DESCRIPTION 


The  schematic  in  Figure  1  illustrates  SHARC’s  module  structure  and  overall 
calculational  sequence.  The  input  module  queries  the  user  for  parameters  needed  to  define 
the  calculation.  Atmospheric  temperature  and  species  density  profiles  are  specified  via  an 
external  file.  Profiles  are  required  for  the  IR-active  species  (NO,  C02,  H20,  03,  CO, 
OH,  and  CH4),  the  major  atmospheric  species  to  which  they  are  collisionally  coupled  (N2, 
02,  O),  and  atomic  hydrogen,  whose  reaction  with  03  provides  the  main  source  of  OH(v). 
Other  input  parameters  include  the  LOS  specifications  and,  if  desired,  the  coordinates  of  a 
localized  auroral  region  through  which  the  LOS  may  pass. 

Excited  vibrational  state  populations  are  calculated  in  the  chemical  kinetics  and 
radiative  transfer  modules  and  are  saved  in  an  output  file  for  later  use.  For  the  quiescent 
atmosphere,  steady-state  conditions  are  assumed.  If  auroral  excitation  is  specified,  a 
time-dependent  chemical  model  calculates  the  additional  production  of  NO  and  NO+ 
arising  from  interactions  of  auroral  electrons.  To  generate  the  desired  LOS  spectrum,  the 
vibrational  state  populations  are  fed  to  the  spectral  radiance  module,  which  outputs 
radiance  and  transmittance  spectra  and  in-band  intensities. 


QUIEUIS. 


PROFILE  OF 

EXCITEO-STATE 

POPULATIONS 


LAYER  PATH  SEGMENTS 


SPECTRAL  RADIANCE 
AND  TRANSMITTANCE 
IN-BAfC  RAOIRNT 
INTENSITY 


Figure  1 .  SHARC  Module  Structure  and  Calculational  Sequence. 
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Condensed  descriptions  of  the  calculations  are  given  below.  Detailed  descriptions 
of  SHARC  may  be  found  in  several  technical  reports.4-6’11,12  Code  upgrades  in  progress 
are  indicated  by  the  dotted  lines  in  Figure  1.  These  upgrades  include  provisions  for  LOS’s 
through  multiple,  distinct  atmospheric  regions,  for  use  in  modeling  terminator,  tidal, 
gravity  wave,  or  other  atmospheric  inhomogeneities. 

2.1  Quiescent  Chemical  Kinetics 

The  quiescent  chemical  kinetics  are  handled  primarily  through  separate  reaction  sets 
for  each  IR-active  molecule;  the  reaction  sets  can  be  solved  selectively  for  the  specific 
molecules  under  investigation.  The  reaction  database  files  contain  the  list  of  vibrational 
states  for  each  molecule,  the  chemical  equations  written  in  symbolic  form,  and  the  rate 
constant  expressions.  The  equations  include  chemical  formation  of  excited  vibrational 
states,  collisional  deactivation  and  excitation  (satisfying  detailed  balance),  spontaneous 
emission,  and  the  radiative  excitation  processes  associated  with  absorption  of  radiation 
from  the  sun  and  from  the  atmosphere.  The  size  of  the  database  is  indicated  in  Table  1 . 
The  atmosphere  is  assumed  to  be  horizontally  uniform  and  is  divided  into  altitude  layers 
that  are  typically  2  km  to  10  km  in  height. 

The  SHARC  CHEMKIN  module,  taken  from  the  Sandia  CHEMKIN 
general-purpose  chemical  kinetics  code,13  reads  the  reaction  files  and  sets  up  the 
time-dependent  differential  rate  equations.  Each  quiescent  vibrational  state  number  density 
[M*]  in  each  atmospheric  layer  is  obtained  from  the  solution  to  the  steady-state  equation 

0  =  d[M*]/dt  =  -A[M*]  +  £sk  -  2>M*]  +  rs+re+ra  ,  (1) 

k  k 

where  A  is  the  total  spontaneous  emission  rate,  the  Sk  are  the  source  terms  for  production 
by  collisional  processes,  and  the  Qk  are  the  collisional  quenching  rates.  rs,  re,  and  ra  are 
the  production  rates  due  to  radiation  from  the  sun,  earth  (i.e.,  from  below  the  50  km 
altitude  boundary  of  SHARC),  and  atmosphere  (above  50  km);  they  depend  implicitly  on 
the  densities  of  the  lower-energy  states  of  species  M.  Equation  (1)  is  solved  for  [M*]  with 
the  constraint  that  the  total  number  density  for  the  molecule  equals  that  specified  in  the 
atmosphere  input  file. 
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Table  1.  Quiescent  Radiating  Species  in  SHARC. 


Molecule 

Isotopea 

States 

Reactions 

Bands 

C02 

626 

28 

217 

34 

co2 

636 

28 

217 

34 

co2 

628 

28 

217 

34 

h2o 

161 

8 

32 

14 

h2o 

181 

8 

32 

8 

h2o 

171 

8 

32 

7 

h2o 

162 

8 

32 

7 

03 

666 

30 

176 

45 

NO 

46 

3 

5 

3 

CO 

26 

3 

9 

3 

OH 

61 

10 

32 

24 

ch4 

21111 

12 

52 

13 

a  1  =  1H,  2  =  2H  or 

12C,  3  =  13C,  4 

=  14N,  6  = 

160,  7  =  170,  8  =  180 

While  most  of  the  excitation  and  de-excitation  processes  involve  only  a  single 
vibrationally  excited  molecule,  leading  to  uncoupled,  linear  equations,  the  species  C02(v), 
H20(v)  and  OH(v)  are  indirectly  coupled  to  each  other  via  resonant  energy  transfer 
processes  involving  N2(v=l),  as  discussed  by  Kumer  and  co-workers.14*16  The 
steady-state  equations  for  these  species  are  linearized  by  equating  the  ground  vibrational 
state  number  densities  with  the  total  number  densities,  and  are  solved  for  the  N2(v=l) 
population.  The  result  is  then  inserted  into  the  uncoupled  Eqs.  (1)  from  which  refined 
excited  state  populations  are  calculated.  For  energy  transfer  processes  involving 
02(v=l),9  02  is  assumed  to  remain  in  LTE. 

The  collisional  rate  constants  were  obtained  from  the  recent  literature  and  from 
Taylor’s  review17  of  measurements  prior  to  1974.  Further  details  are  given  in  a  technical 
report.12  The  solar  excitation  rates  are  derived  from  transmittance  calculations  using 
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either  MODTRAN2  or  the  LOS  radiance  model  discussed  in  Subsection  2.5.  The 
earthshine  rates  kg  are  expressed  in  terms  of  an  effective  blackbody  temperature 
corresponding  to  the  altitude  where  the  vibrational  band  becomes  optically  thick  in  a  nadir 
view. 

The  calculation  of  the  excitation  rate  ra  due  to  atmospheric  emission  originating 
within  the  50-300  km  altitude  range  is  considerably  more  complicated  than  for  rs  and  re, 
and  is  performed  in  a  subroutine  dubbed  NEMESIS,  briefly  described  below. 


2.2  Radiative  Excitation  from  Atmospheric  Layers 


The  effect  of  the  atmospheric  excitation  term  ra  is  to  enhance  the  excited  state 
number  density  within  each  atmospheric  layer  by  an  amount 

=  [M*]-[M*]0  .  (2) 

Here  q  is  the  total  quenching  rate  and  [M*]0  is  the  steady-state  solution  with  ra  set  to  zero, 
as  obtained  from  an  initial  CHEMKIN  run.  Re-excitation  due  to  transfer  from  N2(v  =  l)  is 
incorporated  as  a  reduction  in  q  rather  than  as  a  separate  source  term  in  Equation  (1). 

The  excited-state  enhancement,  and  hence  ra,  may  be  evaluated  following  the 
treatment  of  Kumer  and  co-workers.14,15  For  a  given  vibrational  band  and  atmospheric 
layer  i  we  may  write 

[M*J  =  [M*J0  +  ai«i^Pjiuj[M*j]  (3) 

j 

where  a;  is  the  probability  for  absorption  of  a  photon  entering  the  layer,  P-  is  the 
probability  that  a  photon  emitted  from  layer  j  will  be  absorbed  in  layer  i,  and  gj,  is  the 
branching  ratio  for  re-emission, 


GJj  = 


A* 

A+qi 


(4) 


where  A*  is  the  Einstein  coefficient  for  emission  in  the  given  band. 

For  a  calculated  probability  matrix  (P)  and  set  of  the  set  of  Equations  (3)  for  all 
layers  may  be  solved  simultaneously.  A  physically  instructive  solution  is  provided  by 
repeated  substitution  into  Equation  (3)  of  the  corresponding  equation  for  the  layers  j .  After 
collecting  factors  the  result  is  the  infinite  series 
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(5) 


[M*d  =  [M*j]0  +  aiu,,^(Pji<1>+Pji<2>+...)uj[M*j]0  . 

j 

Here  the  (P)^  matrices  represent  probabilities  for  k-th  order  photon  scattering;  (P)(1)=(P). 
Except  in  the  case  of  the  C02  band,  which  is  strongly  trapped,  the  series  converges  in  a 
small  number  of  terms. 

Both  (P)  and  the  depend  on  the  lower  vibrational  state  populations.  Therefore, 
the  overall  solution  proceeds  in  stages  for  successively  higher-energy  states.  The  first 
CHEMKIN  run  (with  all  ra  =  0)  defines  the  ground  state  populations  [Mj]  and  the  initial 
excited  state  populations  [M*j]0  for  transitions  to  the  ground  state.  The  results  are  used  in 
the  NEMESIS  subroutine  to  evaluate  (P)  and  a,  and  compute  the  radiative  excitation  rates 
ra  for  those  transitions.  CHEMKIN  is  then  rerun  including  the  ra  to  generate  the 
corresponding  populations  for  the  next  set  of  vibrational  transitions.  The 
CHEMKIN/NEMESIS  sequence  is  repeated  until  solutions  for  the  highest-energy  states  are 
obtained. 

The  major  calculational  effort  is  in  computing  the  elements  of  the  layer-layer 
excitation  probability  matrix  (P).  This  involves  a  multidimensional  integral  over  the 
location,  direction,  vibration-rotation  line,  frequency  location  within  the  line,  and 
propagation  distance  of  the  radiation.  The  integral  is  evaluated  with  the  aid  of  Monte 
Carlo  sampling  using  trial  "photons".  The  calculation  assumes  semi-infinite  plane-parallel 
geometry  and  uses  the  exact  temperature-  and  pressure-dependent  Voigt  lineshape.  A  more 
detailed  description  of  the  calculation  is  given  elsewhere.11 

2.3  Sample  Results  for  COp 

Sample  SHARC  results  for  vibrational  states  excited  by  both  atmospheric  and  solar 
radiation  appear  in  Figure  2,  which  displays  vibrational  temperatures  for  selected  states  of 
C02  at  a  solar  zenith  angle  of  82°.  The  states  are  labeled  using  the  HITRAN 
convention;18  only  the  main  isotopes  are  shown.  This  calculation  was  performed  to 
simulate  emissions  observed  near  the  dawn  terminator  during  the  SPIRE  rocket 
experiment. 19  The  species  number  density  and  temperature  profiles  were  obtained  using 
the  MSISE-90  model20  and  an  altitude-dependent  C02  mixing  ratio  representative  of  rocket 
and  satellite  observations.7,21  Up  to  10,000  Monte  Carlo  "photons"  per  band  were  used  in 
the  NEMESIS  runs,  which  was  sufficient  to  essentially  converge  the  vibrational 
temperatures. 
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VIBRATIONAL  TEMPERATURE  (K) 

Figure  2.  Vibrational  Temperature  Profiles  for  Selected  C02  States. 

As  seen  in  Figure  2,  the  ^-excited  states  (00011),  (01111),  and  (10011)  depart 
from  LTE  around  60-70  km  altitude.  Very  similar  vibrational  temperatures  have  been 
reported  by  Wintersteiner  et  al.22  using  the  Air  Force  Phillips  Laboratory  "RAD" 
algorithm.  The  (01101)  (v2)  state  has  been  treated  theoretically  in  recent  studies  by 
Lopez-Puertas  et  al.23  and  Wintersteiner  et  al.7  The  SHARC  calculation  for  Figure  2  used 
the  value  3x10'13v't  cm3/molec/s  for  the  rate  constant  for  collisional  excitation  of 
(01 101)  by  O  atoms,  consistent  with  those  studies,  and  consequently  generates  a  similar 
vibrational  temperature  profile,  which  remains  close  to  LTE  up  to  100  km. 

2.4  Auroral  Kinetics  Model 


The  SHARC  auroral  kinetics  model  calculates  the  enhancements  of  C02,  NO,  and 
NO+  radiation  resulting  from  energy  deposited  in  the  upper  atmosphere  by  solar  electrons. 
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Approximately  600  time-  and  energy-dependent  rate  equations  are  used  to  calculate  the 
secondary  electron  distribution  and  the  subsequent  reactive  and  energy-transfer  processes. 
The  energy  deposition  model  for  the  primary  electrons  is  based  on  work  by  Grun,  Rees  and 
Strickland24"26  as  implemented  in  the  Air  Force  Phillips  Laboratory  code  A  ARC. 27  The 
chemical  reactions  and  energy  transfer  processes  are  described  in  detail  in  References  12 
and  27.  The  integration  of  the  time-dependent  differential  equations  is  accomplished  using 
the  Gear  algorithm  supplied  with  the  Sandia  CHEMKIN  code. 13 

The  use  of  time-dependent  kinetics  results  in  a  significant  improvement  over  the 
steady-state  treatment  used  in  AARC.  For  example,  according  to  SHARC  the  NO(v=l) 
production  efficiency  per  ion  pair  increases  with  time  during  the  aurora  as  the  result  of 
collisional  excitation  of  the  increasing  concentration  of  ground  state  NO. 

2.4. 1  Secondary  Electron  Kinetics 

Secondary  electrons  are  formed  via  collisional  excitation  of  N2,  02,  or  O,  and 
undergo  further  collisions  with  atmospheric  species,  resulting  in  relaxation  of  the  electrons 
and  production  of  ionic  and  electronically  excited  molecular  states.  The  features  and 
altitude  dependence  of  secondary  electron  distributions  in  the  atmosphere  have  been 
discussed  in  detail  elsewhere.28  In  SHARC,  the  secondary  electrons  are  divided  into 
fourteen  energy  bins.  A  distribution  predicted  for  a  Class  III+  aurora  is  shown  in  Figure 
3.  The  strong  dip  around  3  eV  is  due  to  vibrational  excitation  of  N2,  which  in  turn 
provides  the  energy  source  for  the  delayed  C02  v3  auroral  emissions. 

2.4.2  Formation  of  NO  and  NO+ 

When  the  electrons  interact  with  ambient  upper  atmospheric  species,  ionization 
generates  odd  nitrogen  in  the  form  of  N(2P),  N(2D),  and  the  ground  state,  N(4S).  The 
production  mechanisms  in  SHARC  are  direct  dissociation  of  N2  by  the  primary  electrons, 
recombination  of  ions  (e.g.,  NO+)  with  secondary  electrons,  and  charge  transfer  reactions 
of  the  primary  ionic  products.  Collisions  between  nitrogen  atoms  and  02  create  NO.  The 
net  production  rate  is  highly  sensitive  to  the  ratio  of  N(2D)  to  N(4S)  formed,  as  the  NO  is 
formed  primarily  by  reaction  with  N(2D)  and  destroyed  by  reaction  with  N(4S).  The 
vibrationally  excited  NO  is  quenched  by  O  and  02  and  also  undergoes  radiative  relaxation. 
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Figure  3.  Steady-State  Secondary  Electron  Distribution  at  an  Altitude  of  100  km 
for  a  Class  III+  Aurora. 

NO+  is  assumed  to  be  formed  primarily  via  the  ion-molecule  reaction  of  N+  with 
02,  which  competes  with  the  charge  exchange  reaction  forming  02+.  NO+(v)  is  quenched 
by  N2  and  by  radiative  relaxation. 

2.4.3  Vibrational  Excitation  of  C02  and  N2 

CO203)  is  efficiently  excited  by  direct  electron  impact,  which  is  prompt,  and  by 
energy  transfer  from  N2(v=l),  which  is  delayed.  The  vibrationally  excited  N2  is 
presumed  to  arise  from  the  N(4S)  -I-  NO  reaction  and  from  the  quenching  of  0(!D),  and  is 
removed  by  relaxation  with  O  atoms  and  transfer  to  CO203).  In  SHARC,  the  N2(v) 
density  is  calculated  simultaneously  with  the  secondary  electron  distribution. 
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Figure  4.  Prediction  of  the  N2(v)  Production  Efficiency  for  Several  Mechanisms. 

Figure  4  shows  the  N2(v)  excitation  efficiency  for  a  typical  calculation.  It  is  seen 
that  direct  electronic  excitation  of  N2  is  almost  as  important  as  its  chemical  production,  and 
there  is  only  a  slight  altitude  dependence  to  the  efficiencies. 

2.4.4  Radiative  Transfer 

Emissions  from  NO  and  NO+  are  both  optically  thin  and  prompt.  Therefore, 
coupling  of  the  radiative  transfer  and  collisional  processes  is  not  required,  and  the  solutions 
to  the  collisional  mechanisms  suffice.  However,  C0203)  emission  below  100  km  is  both 
optically  thick  and  delayed,  so  that  coupling  of  the  C02  time-dependent  chemical  kinetics 
to  a  radiative  transfer  scheme  is  required.  With  the  present  kinetic  scheme,  a  full  treatment 
would  require  solving  approximately  4000  coupled  differential  equations,  which  is 
impractical.  Therefore,  a  modified  escape  function  approximation  (EFA)  drawing  on  the 
work  of  Kumer29  is  currently  adopted. 
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In  Kumer’s  EFA,  it  is  assumed  that  photons  emitted  from  a  given  layer  either 
escape  the  atmosphere  or  are  re-absorbed  within  that  layer.  A  comparison  between  the 
EFA  and  a  full  NEMESIS  radiative  excitation  solution  for  a  quiescent  atmosphere  is  shown 
in  Figure  5.  Here  "source"  labels  the  C02(v3)  density  calculated  from  CHEMKIN  alone, 
while  the  other  curves  show  the  results  including  radiative  excitation  enhancement.  It  is 
seen  that  the  EFA  agrees  well  with  the  accurate  NEMESIS  calculations  below  70  km  and 
above  130  km;  this  simply  results  from  the  EFA  being  correct  in  the  optically  thin  limit 
and  in  the  limit  of  complete  photon  absorption  within  the  original  layer.  However,  an 
order-of-magnitude  error  can  occur  between  those  altitudes. 

To  insure  that  the  correct  result  is  obtained  in  the  limit  of  negligible  auroral 
excitation,  SHARC  includes  a  correction  which  forces  a  match  to  the  exact  NEMESIS 
solution  for  ambient  conditions.  The  correction  consist  of  a  time-independent  source  or 
sink  term  added  to  the  differential  equation  for  production  of  CO^*^),  whose  magnitude  is 
determined  by  requiring  the  auroral  production  of  C02(f'3)  to  be  zero  at  the  start  of  the 
excitation. 


O 


Figure  5.  Comparison  of  C02(i/3)  Populations  with  Different  Treatments  of 
Atmospheric  Radiative  Excitation. 
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2.5  Line-of-Sight  Spectral  Radiance  Model 


The  LOS  spectral  radiance  is  calculated  by  numerically  integrating  the  radiation 
transport  equation  along  the  LOS  path.  The  approach  in  SHARC  is  based  on  an 
equivalent-width  formulation  for  single  lines.  It  is  around  100  times  faster  than  the  usual 
line-by-line  (LBL)  grid  method,  which  numerically  integrates  over  the  line  shape.  The 
spectral  resolution  of  the  calculation  can  be  made  as  high  as  0.5  cm’1  throughout  the 
50-300  km  altitude  range.  Above  60  km,  the  lines  are  sufficiently  narrow  and  optically 
thin  that  the  resolution  can  be  increased  to  0.1  cm'1. 

The  inhomogeneous  LOS  path  is  represented  with  a  series  of  homogeneous 
atmospheric  layers,  and  the  equivalent  width  is  calculated  using  an  expression  appropriate 
for  Voigt  line  shapes.  The  mean  line  strength  and  width  are  derived  from  the 
Curtis-Godson  approximation.30,31  After  calculating  the  equivalent  width  for  each  line, 
the  combined  equivalent  width  for  all  lines  within  the  spectral  bin  is  calculated  using  an 
approximate  (but  typically  small)  correction  for  line  overlap.  Comparison  to  exact  LBL 
calculations  demonstrate  that  the  overall  errors  in  this  approach  are  less  than  10%  for 
typical  cases. 

2.5. 1  Calculation  for  a  Single  Line 


With  the  atmosphere  represented  as  a  series  of  N  homogeneous  layers,  the  equation 
of  LOS  radiation  transport  for  a  single  molecular  line  is  given  by32 

N 

I  =  Ri  <wi+rw')  -  <6> 

i=  1 

where  the  path  is  from  the  observer  to  the  beginning  of  the  ith  layer,  I  is  the  observed 
radiance,  i  is  the  index  for  each  atmospheric  layer,  R;  is  the  molecular  emission  source 
function  (in  W/cm^sr/cm'1),  and  Wj  is  the  equivalent  width  of  the  line.  The  emission 
source  function  is  given  by 


R 


^l,/3  7 
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where  Cj  is  the  first  radiation  constant  and  y  is  the  ratio  of  upper  to  lower  state 
populations,  defined  by 


Jl_  Pu 

gu  Pi 


(8) 


The  subscript  (1  or  u)  denotes  the  lower  or  upper  state,  g  is  the  degeneracy,  and  p  is  the 
number  density.  The  states  in  Equation  (7)  are  defined  by  their  vibrational,  rotational  and 
electronic  labels,  and  the  number  densities  in  Equation  (8)  are  those  calculated  from  the 
CHEMKIN  and  NEMESIS  routines.  For  LTE  conditions,  y  =  exp(-C2/VT),  where  C2  is 
the  second  radiation  constant,  and  R  reduces  to  the  standard  Planck  blackbody  function. 

The  equivalent  width  of  the  line,  W,  is  defined  as  its  integrated  absorptance.  For  the 
total  N-layer  path,  the  transmittance  within  a  spectral  interval  A v  containing  the  line  is 


r 


1  - 
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_ N 

A  v 
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in  the  limit  that  Av  >  >  WN,  where  WN  is  the  total  path  equivalent  width.  If  Ai> 
represents  the  spectral  bin  width,  then  in  practice  it  need  only  be  slightly  larger  than  WN 
for  the  transmittance  to  remain  physically  valid  (r  <  1)  and  for  the  radiance  I  in  Equation 
(6)  to  be  accurate  for  several-bin- wide  averages.  Within  the  altitude  range  of  SHARC, 
atmospheric  linewidths  are  less  than  0.01  cm'1,  and  equivalent  widths  remain  small  enough 
over  the  entire  spectrum  that  Av  can  be  set  as  low  as  0.5  cm'1. 

For  Voigt  lineshapes,  which  apply  to  upper  atmospheric  conditions,  the 
Rodgers-Williams  approximation33  to  the  single-line  equivalent  width  provides  a 
reasonable  compromise  between  computational  efficiency  and  accuracy.  It  is  given  by 

2  Wi  Wrv 

wv  =  «d^[wE+wJ-(-^)2|  ,  (10) 

where  the  subscripts  V,  D,  L,  and  W  refer  to  Voigt,  Doppler,  Lorentz,  and  weak-line 
limits,  and  aD  (cm'1)  is  the  Doppler  linewidth.  The  equivalent  widths  are  calculated  using 


the  approximations32 
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wjj  =  ln(l+W^,)  , 
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where  a  (cm*1)  is  the  line  width,  S  (cm/molecule)  is  the  line  strength,  and  u 
(molecules/cm2)  is  the  absorber  column  density.  The  errors  in  the  computed  Wv  are 
comparable  to  the  8%  peak  errors  in  the  approximations  for  WD  and  WL.  These 
expressions  for  WD  and  WL  are  more  quickly  calculated  than  the  more  accurate  polynomial 
series  of  Rodgers  and  Williams. 

Since  the  above  expressions  apply  to  a  homogeneous  gas,  and  the  atmosphere 
between  the  observer  and  each  layer  is  inhomogeneous,  Curtis-Godson30,31  path-averaged 
values  of  Su,  aD,  and  aL  are  used.  The  line  parameters  (strengths  and  air-broadened 
half-widths)  are  tabulated  in  a  file  generated  from  the  HITRAN  atlas18  that  has  been 
supplemented  with  lines  for  NO+  and  higher  vibrational  states  of  NO  and  03. 


2.5.2  Calculation  for  Multiple  Lines 


In  the  absence  of  line  overlap,  the  radiances  for  each  line  within  the  spectral  bin 
would  sum  to  form  the  total  radiance  and  the  equivalent  widths  for  each  line  would  sum  to 
form  the  total  equivalent  width.  The  presence  of  overlap  reduces  these  total  quantities.  A 
simple  line  overlap  treatment  is  used  which  provides  reasonable  agreement  with  exact 
results  for  typical  atmospheric  conditions.  The  assumption,  valid  in  the  statistical  limit,  is 
made  that  the  lines  are  randomly  distributed  within  each  bin.  Then  the  combined 
equivalent  width  for  the  total  path,  WN,  is  given  by 


Win  W2N 

WN  =  A,  [  1  -  (1~)(1--^)- "(1- 
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WkN 

Av 


)] 


(14) 


for  "k"  lines  in  the  bin.  Equation  (14)  gives  the  correct  equivalent  width  for  the  r 
calculation  (Equation  (9)).  To  compute  an  approximate  overlap-corrected  radiance,  first 
the  simple  radiance  sum  is  taken  over  all  lines  in  the  bin.  The  result  is  then  multiplied  by 
the  ratio  of  the  combined  equivalent  width  from  Equation  (14)  to  the  equivalent  width  sum 
for  those  lines.  A  more  accurate  formulation  could  be  implemented  if  desired  which  treats 
the  line  overlaps  separately  in  each  path  segment. 
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2.5.3  Illustrative  Calculations 


The  accuracy  of  the  LOS  radiance  algorithm  has  been  explored  through 
comparisons  to  standard  high-resolution  LBL  calculations  for  typical  atmospheric 
conditions.  A  number  of  different  regimes  were  investigated,  including  vertical  and 
horizontal  viewing  geometries,  optically  thin  and  thick  lines,  multiple  overlapping  lines 
and  bands,  and  LTE  and  non-LTE  conditions.  Below,  we  discuss  representative  results  for 
the  C02  15  fim  region. 

Evaluations  of  the  Curtis-Godson  and  Rogers-Williams  approximations  are 
illustrated  by  calculations  on  two  individual  C02  v2  lines,  one  strong  (Q10,  w0  =  667.494 
cm'1,  S  =  3.35xl0'19  cm'Vmolec/cm'2  at  244  K)  and  one  weak  (R56,  co0  =  713.134  cm'1, 
S  =  8.24xl0'22  cm'Vmolec/cm'2  at  244  K).  The  1976  US  Standard  Atmosphere34  was 
used  to  define  the  kinetic  temperature  and  species  number  densities.  The  atmosphere  is 
layered  into  2  km  steps  from  51  to  149  km  and  10  km  steps  from  155  to  295  km.  We 
consider  two  types  of  multi-segmented  paths,  zenith  views  from  60.5  km  altitude  to  65.5, 
71.5,  81.5,  and  101.5  km  altitudes,  and  a  60  km  limb  path.  The  vertical  paths  thus  have 
3,  6,  11,  and  21  segments  along  the  LOS  while  the  limb  path  has  121  segments  along  the 
LOS. 

The  calculated  radiance  values  along  with  their  percent  error  versus  "exact"  LBL 
grid  calculations  are  summarized  in  Tables  2  and  3.  The  errors  range  from  3.4  to  4.9% 
for  the  strong  line  and  from  2.0  to  8.7%  for  the  weak  line.  These  errors  are  consistent 
with  those  expected  for  the  equivalent  width  approximations.  The  limb  calculation  with 
the  strong  line  is  a  good  test  since  the  temperature,  pressure,  and  density  vary  considerably 
along  the  path. 

A  complete  spectrum  for  the  C02  v2  region  is  shown  in  Figure  6  for  a  50  km 
full-limb  view.  This  calculation  provides  a  test  of  the  line  overlap  treatment,  since  at  this 
altitude  there  is  significant  overlap  in  the  Q  branch  of  the  strong  667.380  cm'1  band.  The 
agreement  with  an  "exact"  LBL  calculation,  carried  out  with  FASCOD3,  is  seen  to  be 
excellent  over  virtually  the  entire  300  cm'1  spectral  range  and  the  six-order-of-magnitude 
range  of  radiance.  The  only  major  discrepancies  are  due  to  the  omission  in  this  SHARC 
version  of  three  weak  vibrational  bands,  which  are  responsible  for  the  features  around  508, 
582,  and  757  cm'1. 
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Table  2.  Calculated  Radiances  for  Strong  Line  (Q10)  of  C02  v2. 


FINAL 

RADIANCE 

ERROR 

ALTITUDE 

(lO'8 

W/sr/cm2) 

SHARC 

LBL 

61.5  km 

2.36 

2.34 

0.9 

65.5 

3.25 

3.11 

4.5 

71.5 

3.48 

3.32 

4.7 

81.5 

3.51 

3.35 

4.9 

101.5 

3.51 

3.35 

4.9 

60  km  Limb 

25.9 

26.8 

3.4 

Table  3.  Calculated  Radiances  for  Weak  Line  (R56)  of  C02  v2. 


FINAL 

RADIANCE 

ERROR 

ALTITUDE 

(lO'8  W/sr/cm2) 

_(%) 

SHARC 

LBL 

61.5  km 

0.897 

0.865 

3.7 

65.5 

3.05 

2.81 

8.5 

71.5 

3.54 

3.25 

8.7 

81.5 

3.69 

3.39 

8.7 

101.5 

3.69 

3.39 

8.7 

60  km  Limb 

16.9 

17.3 

2.0 
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3.  DATA  COMPARISONS 


While  comparisons  with  other  computer  codes  can  validate  the  numerical 
procedures,  comparisons  with  upper  atmospheric  IR  data  are  needed  to  test  the  physical 
models  in  SHARC,  including  kinetic  mechanisms,  rate  constants,  and  properties  of  the 
atmosphere.  A  comprehensive  validation  of  these  models  will  require  a  considerable  effort 
and  a  very  large  measurement  database.  Here,  we  present  illustrative  comparisons  of 
predictions  and  data  from  several  rocket  experiments  sponsored  by  the  US  Air 
Force/Phillips  Laboratory. 

Figure  7  shows  a  predicted  emission  spectrum  from  the  quiescent  nighttime 
atmosphere  in  a  near-limb  view  (4°  elevation  angle)  from  an  altitude  of  78  km.  The 
spectral  resolution  has  been  degraded  to  4  cm'1  FWHM.  The  1400  to  3800  cm'1  region 
shown  contains  a  varied  assortment  of  molecular  emission  features,  including  resolved  lines 
of  the  H20  p2  band  (1400-1700  cm"1),  the  NO  Av=l  band  (1700-2000  cm'1),  which  is 
emitted  from  thermospheric  altitudes,  several  03  p3+p3  bands  (2000-2140  cm'1)  generated 
from  chemiluminescent  as  well  as  collisional  processes,  the  CO  Av  =  l  band  (2100-2200 
cm'1),  the  4.3  C02  pj  feature  (2300-2400  cm"1),  consisting  of  a  number  of  optically 

thick,  radiatively  trapped  bands,  and  resolved  OH  Av  =  l  nightglow  lines  (beyond  around 
2500  cm'1).  All  of  these  emission  bands  are  non-LTE. 

The  Figure  7a  calculation  simulates  the  quiescent  atmospheric  radiance  during  the 
recent  EXCEDE  III  experiment,35  which  measured  IR  emissions  induced  by  dosing  the 
atmosphere  with  an  electron  beam.  The  experiment,  flown  in  April  1990  from  White 
Sands  Missile  Range,  NM,  collected  both  beam-on  and  beam-off  data  in  a  slant-viewing 
geometry  using  cryogenically-cooled  Michelson  interferometers  that  covered  a  wide 
spectral  range  (around  500-5000  cm'1)  with  high  resolution  (approximately  2  cm'1).  Late 
in  the  the  flight,  the  sensor  rolled  towards  the  horizon,  and  during  a  beam-off  period  near 
78  km  altitude  it  obtained  data,  shown  in  Figure  7b,  in  an  average  viewing  geometry 
similar  to  that  of  the  SHARC  calculation  in  Figure  7a.  All  of  the  features  predicted  by 
SHARC  appear  to  be  present,  and  they  generally  agree  very  well  in  both  shape  and 
absolute  intensity.  An  exception  is  the  NO  band,  which  is  three  times  stronger  than 
predicted,  indicating  too  low  a  concentration  in  the  SHARC  model  atmosphere.  This 
discrepancy  is  not  very  disturbing  considering  the  large  variability  of  NO.  The  good 
quantitative  prediction  of  the  C02  p3  feature  is  quite  satisfying,  since  it  presents  the  most 
severe  test  of  SHARC’s  atmospheric  excitation  and  LOS  radiation  transport  algorithms. 
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Figure  7.  Comparison  of  SHARC  Calculation  (Top)  and  Observation  (Bottom) 
of  Quiescent  Nighttime  Emission  in  a  Near-Limb  View  from  an 
Altitude  of  78  km. 
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A  more  detailed  examination  of  the  C02  vj  feature  is  provided  in  Figure  8,  which 
shows  limb  radiance  profiles  measured  in  another  experiment,  SPIRE,19  launched  from 
Poker  Flat  near  the  dawn  terminator.  Different  phenomena  dominate  the  excitation  in 
various  altitude  ranges,  and  this  is  reflected  in  the  solar  angle  dependence.  During  the 
day,  the  primary  isotope  is  solar-excited  above  110  km,  while  at  lower  altitudes  the 
radiance  is  enhanced  by  pumping  of  the  2.7  /xm  band  and  by  emissions  from  hot  bands  and 
minor  isotopes.  As  the  density  increases,  the  solar  excitations  are  quenched,  and 
single-quantum  excitation  of  the  primary  isotope  again  dominates.  The  good  agreement 
between  the  data  and  the  SHARC  calculations  indicates  that  these  processes  are  treated 
accurately. 


o 


Figure  8.  C02  4.3  /xm  Band  Limb  Radiance  Measured  in  the  SPIRE  Rocket 

Experiment  and  Predicted  by  SHARC  for  Different  Solar  Zenith 
Angles  Near  the  Dawn  Terminator. 
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Figure  9.  Calculated  and  Observed  NO  Spectrum  for  a  90  km  Zenith  View 
During  a  Class  II  Aurora. 

An  example  of  SH ARC’S  auroral  capability  is  shown  in  Figure  9.  A  nitric  oxide 
spectrum  from  the  1973  Field-Widened  Interferometer  experiment,36  which  observed  an 
IBC  Class  II  aurora,  is  compared  with  a  SHARC  model  calculation  at  a  resolution  of  1 
cm'1.  The  overall  band  shape,  which  includes  strong  aurorally-induced  hot  band 
contributions,  is  reproduced  well  by  the  calculation,  except  near  the  edges  of  the  band 
where  non-LTE  rotational  populations,  not  currently  modeled  in  SHARC,  yield  enhanced 
high-J  lines.  The  overall  absolute  radiance  falls  a  factor  of  four  below  the  data;  however, 
the  results  are  quite  promising  for  an  ab  initio  calculation  that  starts  with  electron  dosing  of 
the  ambient  atmosphere. 

Numerous  other  data/model  comparisons  have  been  and  are  currently  being 
performed  for  all  of  the  major  atmospheric  IR  emission  bands  during  daytime,  nighttime, 
and  auroral  conditions.  The  results  will  be  reported  in  future  papers. 
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4.  CONCLUSIONS 


A  new,  first-principles  computer  model,  SHARC,  has  been  developed  by  the  Air 
Force  for  the  rapid  and  accurate  calculation  of  non-LTE  upper  atmospheric  IR  radiance  and 
transmittance  spectra  with  a  resolution  of  better  than  1  cm*1.  SHARC  treats  the  important 
molecular  bands  from  2  /im  to  40  ftm  (250  cm*1  to  5,000  cm*1)  for  arbitrary  lines  of  sight 
in  the  50-300  km  altitude  range,  accounting  for  the  detailed  production,  loss,  and  energy 
transfer  processes  among  the  vibrational  states.  Calculated  vibrational  temperatures  are 
very  similar  to  results  from  other  non-LTE  codes,  and  SHARC ’s  equivalent- width  spectral 
algorithm  is  found  to  give  very  good  agreement  with  much  more  time-consuming  "exact" 
line-by-line  methods. 

Detailed  comparisons  of  SHARC  radiance  predictions  with  field  measurements  are 
ongoing.  Comparisons  performed  to  date  indicate  satisfying  agreement  for  most  emission 
bands,  including  the  C02  4.3  feature,  which  poses  an  especially  severe  test  of  the 
code.  Data  simulations  using  SHARC  may  be  valuable  in  remote  sensing  of  the  upper 
atmosphere,  such  as  for  deriving  density  profiles  for  variable  species  such  as  NO,  03,  and 
H20. 
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